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Abstract 

If one places N cities randomly on a lattice of size L, we find that Ie^/P and Im-Jv var Y 
with the city concentration p = N/L 2 , where Ie is the average optimal travel distance 
per city in the Euclidean metric and Im is the same in the Manhattan metric. We have 
studied such optimum tours for visiting all the cities using a branch and bound algorithm, 
giving the exact optimized tours for small system sizes (N < 100) and near-optimal tours 
for bigger system sizes (100 < N < 256). Extrapolating the results for iV — > oo, we find 
that Ie\[v = Im\/p = 1 for p = 1, and Ie^/v = 0.73 ± 0.01 and Tm^/p = 0.93 ± 0.02 with 
Im/Ie — 4/7r as p — > 0. Although the problem is trivial for p = 1, for p — > it certainly 
reduces to the standard travelling salesman problem on continuum which is NP- hard. We 
did not observe any irregular behaviour at any intermediate point. The crossover from the 
triviality to the NP- hard problem presumably occurs at p = 1. 
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1 Introduction 

The travelling salesman problem (TSP) is a simple example of a multivariable combinatorial 
optmization problem and perhaps the most famous one. Given a certain set of cities and 
the intercity distance metric, a travelling salesman must find the shortest tour in which he 
visits all the cities and comes back to his starting point. It is a non-deterministic polynomial 
complete (NP- complete) problem [1-3]. In the standard formulation of TSP, we have N 
number of cities distributed randomly on a continuum plane and we determine the average 
optimal travel distance per city Te in the Euclidean metric (with Ave = \^Ax 2 + Ay' 2 ), or T M 
in the Manhattan metric (with Atm = | Aa;| + | Ay|). Since the average distance per city scales 
(for fixed area) with the number of cities N as 1/ v^V, we find that the normalized travel 
distance per city = IeVN or VLm = TmVN become the optimized constants and their 
values depend on the method used to optimize the travel distance. Extending the analytic 
estimates of the average nearest neighbour distances, in particular within a strip and varying 
the width of the strip to extremize (single parameter optimization approximation), one gets 
| < VL E < 0.92 [4] and ^ < VL M < 117 [5]. Careful (scaling, etc.) analysis of the numerical 
results obtained so far indicates that VL E ^ 0.72 [6]. 

Similar to many of the statistical physics problems redefined on the lattices, e.g., the 
statistics of self-avoiding walks on lattices (for investigating the linear polymer conforma- 
tional statistics), the TSP can also be defined on randomly dilute lattices. The (percolation) 
cluster statistics of such dilute lattices is now extensively studied [7]. The salesman's op- 
timized path on a dilute lattice is necessarily a self-avoiding one; for optimized tour the 
salesman cannot afford to visit any city more than once and obviously it is one where the 
path is non-intersecting. The statistics of self-avoiding walks on dilute lattices has also been 
studied quite a bit (see e.g., [8]). However, this knowledge is not sufficient to understand 
the TSP on similar lattices. The TSP on dilute lattices is a very intriguing one, but has not 
been studied intensively so far. 

The lattice version of the TSP was first studied by Chakrabarti [9] . In the lattice version 
of the TSP, the N cities are represented by randomly occupied lattice sites of a two- di- 
mensional square lattice (L x L), the fraction of sites occupied being p (= N/L 2 , the lattice 
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occupation concentration). One must then find the shortest tour in which the salesman 
visits each city only once and comes back to its starting point. The average optimal travel 
distance in the Euclidean metric Ie, and in the Manhattan metric Im, are functions of the 
lattice occupation concentration p [10]. We intend to study here the variation of the nor- 
malised travel distance per city, Qe = Ie-^/P and = ^Mi/P, with the lattice concentration 
p for different system sizes. It is obvious that at p — 1, all the self- avoiding walks passing 
through all the occupied sites will satisfy the requirements of TSP and Q E — 1 — Vt M (the 
distance between the neighbouring cities is equal to the unit lattice constant and the path 
between neighbouring sites makes discrete angles of of n/2 or its multiples with the Carte- 
sian axes). The problem becomes nontrivial as p decreases from unity: isolated occupied 
cities and branching configurations of occupied cities are found here with finite probabili- 
ties and self-avoiding walks through all the occupied cities, and only through the occupied 
cities, become impossible. As p decreases from unity, the discreteness of the distance of the 
path connecting the two cities and of the angle which the path makes with the Cartesian 
axes, tend to disappear. The problem reduces to the standard TSP on the continuum in the 
p — > limit when all the continuous sets of distances and angles become possible. We study 
here the TSP on dilute lattice employing a computer algorithm which gives the exact opti- 
mized tours for small system sizes (N < 100) and near-optimal tours for bigger system sizes 
(100 < N < 256). Our study indeed indicates that Qe and Qm vary with p and Qe — 0.73 
and Vt M — 0.93 as p — > 0. 
2 Computer Simulation and Results 

We generate the randomly diluted lattice configurations following the standard Monte Carlo 
procedure for different system sizes. For each system size N, we vary the lattice size L so 
that the lattice concentration p varies. For each such lattice configuration, the optimum 
tour with open boundary conditions, is obtained with the help of the GNU tsp_ solve [11] 
developed using a branch and bound algorithm (see Fig. 1). It claims to give exact results 
for N < 100 and near-optimal solutions for 100 < N < 256. It may be noted that the 
program works essentially with the Euclidean distance. However there exists a geometric 
relationship between the Euclidean distance and the Manhattan distance. We may write 
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Ie = Yn=i r ii and Im = Hi=i r i a ii where is the magnitude of the Euclidean path vector 
between two neighbouring cities and r jcc, = r\ ( | sin ^ | + 1 cos 9i\) is the sum of the components 
of the Euclidean path projected along the Cartesian axes. Naturally, 1 < on < y/2. If l E 
corresponds to the shortest Euclidean path, then YliLi K > Y^iLi r % > f° r any other path 
denoted by the primed set. If the optimized Euclidean path does not correspond to the 
optimized Manhattan path, then one will have Y^=\ r i a i < Sili r %d%-, where all the ctj 
and a[ satisfy the previous bounds. Additionally, for random orientation of the Euclidean 
distance with respect to the Cartesian axes, (aij ) = (a/) = (2/ti) £ 12 (sm 6 + cos 6)d6 = 4/vr. 
It seems, with all these constraints on a's and a n s, it would be impossible to satisfy the 
above inequalities on J2 r i, an d J2 r i a i- 111 fact, we checked for a set of 50 random optimized 
Euclidean tours for small N (< 10), obtained using the algorithm, whether the optimized 
Manhattan tours correspond to different sequence (of visiting the cities), and did not find 
any. We believe that the optimized Euclidean tour necessarily corresponds to the optimized 
Manhattan tour. We then calculate Ie and Im for each such optimized tour. 

At each lattice concentration p, we take about 100 lattice configurations (about 150 
configurations at some special points nearp — > 0) and then obtain the averages Ie and Im- We 
then determine Qe — Ie y/p and Qm = Im \/p and study the variations of Qe and Qm, and 
of the ratio Qm/^e with p. We find that Qe and Qm both have variations starting from the 
exact result of unity for p — 1 to the respective constants in the p — > limit. In fact we noted 
that although Qm continuously decreases as p — > 0, it remains close to unity for all values of 
p. We studied the numerical results for N = 64, 81, 100, 121, 144, 169, 196, 225 and 256. 
The results for iV = 64 and 100 have been shown in Figs. 2 and 3 respectively. We have 
studied the variations in the values of Qe and Qm against 1/N for p — > 0, to extrapolate its 
value in the iV — > oo limit. It appears that for the large N limit (see Fig. 4), Qe{p — > 0) and 
Qm(p — > 0) eventually extrapolate to 0.73 ± 0.01 (as in continuum TSP) and to 0.93 ± 0.02, 
respectively. This result for Qe (at p — > 0) compares very well with the previous estimates 
[6]. As p changes from 1 to 0, the ratio Qm/Qe changes continuously from 1 to about 
1.27 (~ 4/ it) (see Fig. 4), which is the average ratio of the Manhattan distance between two 
random points in a plane and the Euclidean distance between them [10, 5]. 
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3 Conclusions 

We note that the TSP on randomly diluted lattice is certainly a trivial problem when p — 1 
(lattice limit) as it reduces to the one-dimensional TSP (the connections in the optimal tour 
are between the nearest neighbours along the lattice). Here Qe(p) — &>m{p) = 1- However, 
it is certainly NP- hard at the p — > (continuum) limit, where Qe — 0.73 and Q M ~ 0.93 
(extrapolated for large system sizes N). We note that VL M remains practically close to unity 
for all values of p < 1. Our numerical results also suggest that VL m /&e — 4/V as p — * 0. It is 
clear that the problem crosses from triviality (for p — 1) to the NP- hard problem (for p — > 0) 
at a certain value of p. We did not find any irregularity in the variation of Q at any p. A 
naive expectation might be that around the percolation point, beyond which the marginally 
connected lattice spanning path is snapped off [7], the or VLm suffers some irregularity. 
The absence of any such irregularity can also be justified easily: the travelling salesman has 
to visit all the occupied lattice sites (cities), not necessarily those on the spanning cluster. 
Also, the TSP on dilute lattices has got to accomodate the same kind of frustration as the 
(compact) self-avoiding chains on dilute (percolating) lattices, although there the (collapsed) 
polymer is confined only to the spanning cluster. This indicates that the transition occurs 
either at p = 1_ or at p = 0+. From the consideration of frustration for the TSP even at 
p — 1_, it is almost certain that the transition occurs at p — 1. However, this point requires 
further investigations. 

Acknowledgement : We are grateful to O. C. Martin and A. Percus for very useful com- 
ments and suggestions. 
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Figure captions 



Fig. 1 : A typical TSP for (N =) 64 cities on a dilute lattice of size L = 30. The cities are 
represented by black dots which are randomly occupied sites of the lattice with concentration 
p = N/L 2 ~ 0.07. The optimized Euclidean path is indicated. 

Fig. 2 : Plot of Q E , Vl M and £l M /Q E against p for iV = 64 cities, obtained using the 
optimization programs (exact). The error bars are due to configurational fluctuations. The 
extrapolated values of He, and Qm/^e are indicated by horizontal arrows on the y-axis. 
Fig. 3 : Plot of Q E , tt M and Vl M /^E against p for iV = 100 cities, obtained using the 
optimization programs (exact). The error bars are due to configurational fluctuations. The 
extrapolated values of Q E , and Qm/^e are indicated by horizontal arrows on the y-axis. 
Fig. 4 : Plots of Q E , Qm and of Qm/^e m the p — > limit, against 1/N. The error bars 
are due to configurational fluctuations. The extrapolated value of fl E , Om and VLm/^e in 
this p — > limit for N — > oo are indicated by horizontal arrows on the y-axis. 
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